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I. INTRODUCTION 



One of the main goals of the future experiments at hadron coUiders will be undoubtedly 
the search for supersymmetric particles. In the specific theoretical framework of the 
MSSM, a vast amount of hterature already exists showing the expected experimental 
reaches for the different supersymmetric spectroscopics and parameters, both for the 
TEVATRON [1] and for the LHC [2] cases. For what concerns the possibility of direct 
production, nothing has to be added (in our opinion) to the existing studies, leading to 
the conclusion that, if any supersymmetric particle exists with a not unfairly large mass, 
it will not escape direct detection and identification. 

A nice and special feature of the present and future hadron colliders is the fact that, 
on top of direct SUSY production, a complementary precision test of the involved model 
(assuming a preliminary discovery) is also in principle possible. This implies the study 
of virtual SUSY effects in the production of suitable (not necessarily supersymmetric) 
final states, in full analogy with the previous memorable analyses performed to test the 
Standard Model at LEPl, SLC and LEP2 [3]. For what concerns the realistic experimental 
accuracy requested by a similar search, one expects [4], [5] a possible relative few percent 
level. In this spirit, the existence of virtual SUSY effects as large as a relative ten percent 
or more should be carefully examined and investigated. 

The aim of this preliminary paper is actually that of showing that the production 
of a final heavy (6, t) quark pair at LHC could be particularly convenient for the search 
of virtual SUSY effects. This is due to a number of technical features (cancellations 
of disturbing contributions at high energies, not negligible final quark masses,...) that 
will be thoroughly illustrated in the following Sections. In this preliminary analysis, we 
studied the electroweak and the SUSY QCD contributions to the subprocess initiated by 
qq and gg partons and we limited our application to the invariant mass distribution of 
PP qq + ... at LHC using standard quark and gluon structure functions. The rewarding 
result was that of finding that the SUSY electroweak and the SUSY QCD terms concur 
to produce an overall effect that can be as large as a relative twenty percent in the 
assumed hght SUSY scenario. This effect is more due to the electroweak component than 
to the QCD component, essentially owing to the final heavy quark Yukawa contribution 
(proportional to the quark squared mass). The picture at one loop is in fact, surprisingly 
to us, practically the same that one would find in the case of electron-positron annihilation 
at LC [6], as we shall briefly discuss in this paper for sake of comparison. At the possible 
few percent experimental accuracy level [5], this effect should not escape a dedicated 
measurement, and could provide an important test of the (assumedly -and hopefully- 
discovered) supersymmetric scenario. 

Technically speaking, this paper is organized as follows: Section 2 will be devoted 
to a description of the one-loop analysis for an initial qq state; in Section 3, the same 
description will be given for an initial gg state and a brief comment will be given on the 
analogy of the results with those obtained for an initial electron-positron state; in Section 
4, a numerical analysis of the possibly visible effects on the realistic proton-proton initiated 
process will be shown and a few conclusions will finally appear in Section 5. 
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II. INITIAL qq STATE 



We begin our analysis at the partonic level considering as initial states the light quark 
components of the proton (including the bottom component but ignoring the top quark 
one). At lowest ag order the scattering amplitude is given by the Born terms schemati- 
cally represented in Fig.l. The s-channel diagram of Fig. la applies to all initial qq pair 
annihilation and leads to: 

= EC"'^ nv{<l)^P^u{q)M<l')^PA<l')] (2.1) 

whereas the t-channel diagram will only apply to bb — > bb scattering (because we shall 
neglect the tt — > tt contribution) and writes: 

where i,j refer to or L quark chiralities with Pr l — I is the intermediate gluon 
colour state and 

„Born s _ ^'^'^s „Born t _ '^'^'^s /r, o\ 

O'ij - a^j - K^-'i) 

At the Born level, electroweak contributions only consist in the Drell-Yan process 
QQ ^ 1: Z ^ q'q' but the effect in the cross section will be reduced by a factor a^/ ~ 
0.01 as compared to the QCD one and will be neglected at the expected LHC accuracy. 

Our starting assumption, as already stated in the Introduction, will be the discovery 
of (at least some) supersymmetric particles. In this spirit, we shall examine the possibility 
of performing a precision test of the candidate model by looking at its higher order virtual 
effects in the considered heavy quark production at hadron colliders, starting from the 
determination of these effects at the partonic level. As a first case to be examined, we 
shall restrict our study to the MSSM framework. For what concerns the higher order 
diagrams to be retained, we must at this point make our strategy completely clear. In 
general, the higher order corrections to the considered processes will be of electroweak and 
of strong nature. We shall begin our analysis with the detailed study of the electroweak 
components given in the following Section Ila. 



A. Electroweak corrections 

To the first order one-loop level the electroweak corrections will typically correspond 
to the diagrams schematically represented in Fig. 2. Note that we considered in the elec- 
troweak sector both the SM components and the genuinely supersymmetric ones, that 
will be often denoted for simplicity with a (SUSY) apex. In our analysis, wc shall not 
consider higher order (e.g. two loops) electroweak corrections. This (pragmatic) attitude 
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is, at least, supported by the fact (that we showed in a previous paper [7]) that such terms 
at lepton colhders (LC) are only requested if the available cm. energy of the process is 
beyond, roughly, the 1 TeV range. In our investigation of proton colliders, the structure 
of the leading electroweak corrections is essentially similar to the LC one (with straight- 
forward modifications of the initial vertex contributions) . To reproduce the benefit of the 
absence of hard resummation computations, we shall therefore be limited to cm. energy 
values not beyond the 1 TcV limit. 

Having fixed the considered energy range, a numerical evaluation at one loop can 
now be performed in a (reasonably) straightforward way. A great simplification can 
nonetheless be achieved under the assumption that the cm. energy is sufficiently larger 
than all the masses of the (real and virtual) particles involved in the process. In this case, 
a logarithmic expansion of the so-called Sudakov kind can be adopted. All the details of 
such an approach for the specific MSSM model have been already exhaustively discussed 
in previous references for an initial electron-positron state [8], [9], but the results are 
essentially identical if the initial state is a qq pair (as we said, one must only modify 
in a straightforward way the initial vertex contribution). To avoid loss of space and 
time, we defer thus the reader to [8], [9] for details. The only relevant point that wc 
want to make is the fact that the "canonical" logarithmic expansion is given to next-to 
leading order, i.e. retaining the double and the linear logarithmic terms, and ignoring 
possible extra (e.g. constant) contributions. To this level, the remarkable simplification 
for the considered heavy quark final state is that only two SUSY parameters appear in 
the expansion. One parameter is tan/5, contained in the coefficient of the linear In of 
Yukawa origin; the second one is a (common) SUSY mass scale M, to be understood as 
an "average" supersymmetric mass involved in the process. 

Given the fact that the supersymmetric masses might be not particularly small (the 
existing limits for squark and gluino masses arc at the moment consistent with a lower 
bound of approximately 300 GeV [1], [10]) and taking into account the qualitative 1 TeV 
upper bound on the cm. energy requested for the validity of a one-loop expansion, we 
shall adopt in this preliminary analysis a reasonable working compromise. In other words, 
we shall consider energy values in the 1 TeV range, and a SUSY scenario "reasonably" 
light i.e. one in which all the relevant SUSY masses of the process are smaller than, 
approximately, 350 — 400 GeV. This will encourage us to trust a conventional Sudakov 
expansion, supported by previous detailed numerical analyses given for electron-positron 
accelerators [11]. 

The previous discussion was related to the SUSY electroweak effects. Our attitude 
will be in conclusion that of considering both the SM and the SUSY components at the 
same one loop level. Strictly speaking, since we are only interested in the SUSY effect, 
we could even ignore the SM contribution. The reason why we retain it is simply that 
it is easy to estimate it and interesting to show the "SUSY enhancement" in the MSSM 
overall correction. 

After this long but, we hope, sufficiently clear discussion, we are now ready to illus- 
trate the one-loop numerical details. With this aim, we shall start from the electroweak 
diagrams of Fig. 2 and write, following our previous notations [7-9]: 



4 



Electroweak logeirithmic corrections on the amplitudes 



For each qq q'q' amplitude we can write, following ref. 

j^l loop ew A^'^^^^\ -\- C^^] 

with 



(2.4) 



C — Cjn, gauge ~l~ Qn, Yuk ~l~ ^fin gauge ~l~ ^fin, Yuk ~l~ Q 



ang 



(2.5) 



where Ci„ refers to the universal (i.e. process independent) corrections due to the initial qq 
pair, Cfin to the universal corrections due to the final bb or tt pair and Cang is a process and 
angular dependent correction typical of each qq q'q' case. The Yukawa contributions 
are sizable only for q (or q') being 6 or t quarks. The coefficients c'^^ receive SM and SUSY 
contributions, depending on the quark chirality. The SM contributions are {u represents 
both up and charmed quarks, and d both down and strange quarks): 



„uu SM 
"in, gauge L 

„uu SM 
"in, gauge R 



tt SM 

fin, gauge L 



tt SM 

fin, gauge R 



a{l + 26c^^) 



[3/n 



Mir 



-In' 



a 



Mir 



-In' 



Mir' 



(2.6) 



^dd SM 
"in, gauge L 

^dd SM 
"in, gauge R 



^bb SM 
"fin, gauge L 



C 



bb SM 

fin, gauge R 



a(l + 26c^) 
1447rs^c^ 

[3/n 



a 



[3ln 



367rc^ 



Mir 



Mir 

-In' 



-In' 



Mir' 



Mir' 



(2.7) 



J)b SM 

"in, fin, Yuk L 



— a 



ft SM 

in, fin, Yuk L 



J)b SM 

"in, fin, Yuk R 



Jt SM 

"in, fin, Yuk R 



a 



a 



a 



w 



2 



][ln 
][ln 



Mir' 



Mir' 



Mi; 



+ 



Mir 



][ln 



M2 ^ 



(2.8) 



This last angular dependent term refers to chirality states {ij) = {LL, LR, RL, RR) 
and involves the quark charges Qg and the Zqq couplings g^j^ = /gi(2 — 4|(5q|s^), 

9qR = ^'^QqSw- 
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The additional SUSY contributions affect only the universal parts: 

uu SUSY _ tt SUSY _ _ Q^(l + 26c^) s 

^in, gauge L ^fin, gauge L 14:4:71 S^^C^y 

uu SUSY _ tt SUSY " [/„ 1 (2 10) 

'-in, gauge R '-fin, gauge R dnC^^ M'^ y"^'^^) 



dd SUSY ^ bb SUSY ^ _ Q:(1 + 26c^) S 

'-in, gauge L '-fin, gauge L 14471 Sy^Cy^ 

dd SUSY ^ bb SUSY ^ _ ^ \in ^ ] (')^^\ 

Hn, gauge R fin, gauge R SGtTC^ \'^-'-^) 



a r m 



2 



<?f.f^y.. R = - 8^[]i^(l + 2cot^/3)][/n^] (2.12) 



One can see (as emphasized in ref. [9]) that the total MSSM gauge correction is 

quite simply obtained by replacing the logarithmic SM factor [Slrijp ^^^a^] 

[2lnjp In^jp-] and the total MSSM Yukawa correction by replacing the SM by 

2m(2(l + cot^ p) and ml by 2mg(l + tan^ /5). 



Results for angulcir distributions, averaged over initial, summed over final 
polarizations 

We now list exphcitely the complete MSSM results for each subprocess. According to 
the rules given just above it is easy to separate the pure SM and the additional SUSY 
parts. We first consider the subprocesses involving only the annihilation channel of Fig. la 
and Fig.2, and write: 

J 1 loo-p fjfj^orn 2 

=^^^^ + !l!f±{ (l + cos'm [5]+2cos^ [D] } (2.13) 
dcose dcosO 18s ^ ^ ' ^ ^ ^ ^ ^ ^ ' 

with 



7 ;: — -^r^ ( 1 + cos 

and 



dcos9 9s ' ^ ^ ^ 
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For uu — > bb 



For dd — > bb 



For uu — > tt 



For (id — > tt 



'ddtt 



-£|;[|;(l + eot^/3)H-g(l + tar.^/.)]R.^] 

- ^i%^/ni^^[/n^] (2.15) 
IStts^c^ 1 + cos^^ M2,J ^ ^ 

a 1 — cos ^ ,. s , 



Ot 



q;(8sL — 9), 1 — cos^., s , 

+ 2 2 Ti^^TT^ 2.17 

IStts^c^ 1 + cos^^ M^^ ^ ^ 



Sum = 2 2 (27 - 10s^)[2/n— ^ - /n^ 



187rs|^.cfj. 1 + cos^^ M^^ ^ ' 



^ a , 1— cos^r, s , 

27rs^c^ 1 + cos e 



Gi S S 

3^(27- 16.f,)[*^-ta^^l 



a |.3mj 2 



2 



(1 + cot' /?) + TTF (1 + tan' /?)] [in- 



47rs^^M^' My ''^"^ My 

„(4.5,-9),^l-oos.|^ . J (2.21) 



l%T:syi^ l + cos^^ 
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„ a , 1 — cos^ s , , , 



For bb tt 



- T^[S^(l + cot'/5) + ^(l + tan2/3)][/n-^, , 
At^s^M^^ My ^^^^ 

a(44.-9) 1-cosg g 



^ a , 1 — cos6' s , ^ ,s 

A6tt= ^ 9 9 fa:; niln—pr] (2-24) 

In the si)e( ial case of bb hb we have to add the annihilation amphtude of Fig. la and 
of Fig. 2, and the scattering amphtude of Fig.lb and of the s — > i crossed diagrams of 
Fig. 2. This leads to the result 



da^ ^""Pibb bb) da^^^^'ibb bb) 



a 



w 



' ' " ' cot^/?) + — 2^(1 + tan2/3)][/n- 



ct^a rj s , , , 18 — 8s^, , i-M^ , 1 — cos^ 

" T^^^'^A^J ^ ^ ds^^cl, ~ 3^^ 1 + cos^^ 

u"^ 1 + COS6' 27 - 22s^ _^ ti^ l-cos6' 

37t^ 2 > + 9s2,c2^ 3^^''' 2 

, 4.^ ^,,1^ _ ^ ,,1±^ ) y (2.25) 



with 



9sf^.c|^.'' s2 i + cos^ t^ 2 



d(7^°^"(6& ^ 66) 27ra2 ^ t2 ^2 ^ ^2 21^2 



dcos9 9s ' ^ ^2 35 J (2-26) 

In the expression (2.25), the first part (which factorizes the Born term) is the universal 
effect including the gauge term and the double Yukawa effect, whereas the second part 
is the angular dependent effect from s and from t channels (one can check, restricting to 
the terms, that one recovers the previous dd ^ bb case); the s — > t crossing relation 
between annihilation and scattering contributions is also clearly satisfied. 

Having completed the discussion of the electroweak SUSY effect, we now move in the 
following subsection lib to the discussion of the strong (QCD) SUSY correction. 
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B. QCD SUSY correction 



A preliminary statement to be made before entering the discussion of the QCD SUSY 
corrections is that we shall treat this effect under the same assumptions that we adopted 
for the electroweak sector, subsection Ila. In other words, we shall still concentrate our 
analysis on cm. energy values in the 1 TeV region, assuming the previously considered 
"reasonably" light SUSY scenario. This will allow us to use the same kind of simple 
logarithmic expansion that was exploited in the electroweak case, at the one loop level. 

For what concerns the expected validity of a one-loop perturbative expansion, the 
situation is now, though, drastically different from that of subsection Ila, and requires a 
precise subtle choice of strategy. It is actually well-known [12] that, in the SM domain, 
the one-loop QCD correction is not accurate enough and higher order terms seem to be 
fundamental. This is understandable in an extremely simplified fashion as a consequence 
of the small scale which enters the SM running of ctg. In this qualitative picture, one 
would expect that the corresponding SUSY effects do not share this dramatic problem, 
owing to the much larger mass scale involved. This would justify the expectation that, for 
the restricted subset of SUSY QCD corrections, a one-loop calculation can be sufficient. 

Another essential difference between the SM QCD and the SUSY QCD corrections is 
that in the SM part the infrared logarithms arising from virtual gluon contributions cancel 
against those occurring in soft real gluon emission, leaving only "constant" terms. In the 
SUSY case, large logarithms of virtual origin and scaled by the average SUSY mass will 
remain. A practical attitude seems to be therefore that of isolating the SM higher order 
QCD effects, considering them as a "known" quantity, much in analogy with what is done 
for the canonical QED corrections, and to factorize an explicit one-loop term containing 
the genuine SUSY QCD correction, to be added in the usual one-loop philosophy to the 
electroweak one. This is the approach that we shall follow in the rest of the paper, which 
is, as we said, only interested in the evaluation of the genuine overall SUSY effect. 

Having made this preliminary statement, we move now to the evaluation of the QCD 
SUSY effects at one loop. These are represented schematically in Fig. 3, and can be 
classified in two quite different sectors, respectively of vertex kind and of RG origin. The 
first ones are essentially similar to analogous electroweak vertices of Sudakov kind, with 
a simple replacement of gauginos by gluinos. They produce, in the adopted "Sudakov 
regime" , the following effect on the amplitude for each external qq pair 

,,SUSY QCD ^ _ S 

which leads to an additional coefficient in the series eq. (2.4,2.5) equal to twice this value 
for Sk qq ^ q'q' annihilation amplitude or for a, qq ^ qq scattering amplitude, i.e. eq.(2.4) 
is replaced by 

^Born[i^^ew_ ^^n^] (2.28) 

From eq.(2.28) a rather important (in our opinion) feature can be stressed. Numerically, 
one sees that for "reasonable" values of the squark and gluino masses, the QCD SUSY 
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vertex effect has a numerical value of approximately 5 percent at 1 TeV cm. energy. 
This supports our assumption that higher order terms can be neglected, as we shall do 
in this paper. Note also that the size of the effect is smaller than that of the analogous 
electroweak component, particularly for large tan/3 values in the Yukawa couplings. In 
other words, at the one loop level, electroweak SUSY corrections appear to be "stronger" 
than the corresponding QCD ones, in full analogy with a similar feature first stressed in 
the electron-positron case [13]. 

From a practical point of view, the very welcome feature that characterizes the QCD 
SUSY vertex is the fact that the sign of the one-loop effect is the same (negative) as 
that of the corresponding electroweak one. As a consequence of this rewarding concur- 
rence, the overall SUSY vertex effect at one loop gets an enhancement that will improve 
the possibility of experimental detection, and we shall return to this point in the final 
discussion. 

The next term to be computed corresponds to the RG diagram of Fig.(3e). Following 
our previous discussion, we shall only consider the genuine SUSY effect on the intermediate 
gluon bubble. This is given at one loop from standard formulae [12]. In a super-simplified 
attitude of considering the relevant cm. energy beyond all SUSY scales, it would lead to 
the following modification: 

Saf'^is) = a,iM^)[-Bsusy^^ln^,] (2.29) 

The SUSY contribution arising in the MSSM for s > is obtained using BsusY — 
—2. This leads effectively to one more additional coefficient in the series of corrections 
eq. (2.4,2.5) to the annihilation amplitude 

qq,SUSYQCD „ tts(^^), S (n 0(\\ 

Crg = -BsusY /n— (2.30) 

and a similar correction to the scattering amplitude with ln{s/M'^) replaced by 
ln{-t/M^). 

Note that, in a less optimistic situation of larger SUSY masses, the effect would be 
reduced and should be computed more carefully. But, independently of this, we notice 
a less pleasant feature of eq.(2.30) compared to eq.(2.28) i.e. the fact that it produces 
an effect of opposite sign with respect to that of the electroweak one, thus reducing the 
overall SUSY correction: 

j^l loop SUSY _ ^fiorn j-^SUSY ew _|_ ^SUSY QCDj ^-^^ 

with 

^susY QCD ^ _ '^inJ- + ^/n4- = ^/n-^ (2.32) 

37r M2 TT M2 37r M2 ^ ^ 

From this point of view, the qq initial state exhibits a "disturbing" feature for the 
detection of SUSY effects. This feature will not affect the determination of the SUSY 
correction for the initial gluon gluon state at high cm. energies, as we shall show in the 
following Section 3 which will be devoted to the study of that process. 
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III. INITIAL gg STATE 



The amplitude for the process g'^g^ — > q'q' (where i,j denote the gluon color states) is 
obtained by summing the 2 diagrams of Fig.4a,b and the crossed diagram of Fig.4b: 

Af"^- = -i ^^r\e\e^)[u{q')^nk' - y)An (3-1) 
A^"'-" = - f [^(?')^(7'^4)(7''(^^ -/)p)(7^4)^(g")] (3.2) 

^Born ^ _ ^[^(g')--(7^6^j(7''(fc^ -p^"')p)(7^6L)^(g')] (3-3) 

where s ^ {k' + = +p«')^ t = {k'-p^'Y, u = {k'-p^'f, and {e\ k'), (e^ k^) are 
the polarization vectors and four-momenta of the gluons. 

It is important to notice the gauge cancellations occurring at high energies. From the 
above expressions one can immediately compute the helicity amplitudes 
with r* = ±l, r'' = ±l, A"^ =±1/2, A"^ =±1/2 being the gluons and quark helicities. 
At high energy, neglecting quark masses, one obtains only chirality conserving terms with 
A^' = — A^' = A = ±1/2. The contribution to the amplitudes with = —r^ = r = ±1 
arises from t and u channel terms: 

mBornr x 2 / '^^'^^ x 2A COS ^ ± T . _ , 2/ -^^-^N 2A COS ^ ± T . „ 

F«-(.,-.,A,-A) =9f( — )^-^.m9 + 9!(— )^-^^.m« (3-4) 

whereas the = r-' = r = ±1 amplitudes get contributions from s, t and u channel 
terms: 

F^-(r,r,A,-A) = -z^,V^^-^^(2A) sin^ ± ^,2(::^)(2A) sin^ - ^2^^^(2A) sin^ (3.5) 

and totally cancel as — = i/*-'^^. 

So, at high energy, we are only left with contributions to F^'"""(r, — r. A, —A) arising 
from the t and u channel quark exchange diagrams, A = ±1/2 corresponding to R,L 
chiralities respectively. 

The electroweak corrections are then extremely simple. At first order no electroweak 
correction arises for gluons. Only the universal gauge and Yukawa terms appear for the 
final b or t quark pair: cf£^ and 4^^ ^^^ from eq.(2.5). 

The SUSY QCD corrections turn out to be also extremely simple. Because of the 
gauge cancellation of the s channel term, at first order we can ignore the SUSY QCD 
corrections to this part and only consider the t and u channel quark exchange diagrams. 
Fig. lb and Fig. 5. There is no SUSY QCD correction to the external gluon lines because 
of the cancellation between the gluon splitting function and the gluon coupling Parameter 
Renormalization (this fact is similar to the one occurring for electroweak gauge bosons as 
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noticed in ref. [14,15]). Only the universal SUSY QCD correction to the external quark 
lines appear, given by c[^\^^^^ '^'"^ of eq.(2.27). 

The angular distributions averaged over initial gg, summed over final b, b or t,i polar- 
izations are then given by 

da^ ^""P(gg ^bb) da^°'''{gg ^ bb) _ « 2 ^r^, « , 2 « n 

- { 1 + T77-^(27 - 22.^)[2Zrz-^ - In' ' 



dcosO dcosO ' U47lsl^,ci^ " M^' 

" [^il + com + '<il + tanmin^]-'^[ln4,]} (3-6) 



da^ ^°"P{gg ti) da^""'{gg tt) 



a 



{ 1 + 7773^(27 - 10.^)[2/n-^ - In' 



dcosO dcosO ' 1447rs2^,c^^ 

- ^[^(^^^m^^il^tanmin^\- '-^[ln^\ } (3.7) 

In the above expression we have written the complete MSSM electroweak correction 
and the SUSY QCD correction. In the MSSM electroweak part one can easily separate 
the pure SM and the SUSY components using the rules already stated in the previous 
section. In particular, the SUSY modification of the universal terms consists of replacing 
the linear SM ~ 3/n term by 2/n, thus leading to a negative —In effect, while the Yukawa 
term produces a substantial negative contribution due to the tan j3 parameter. 

The Born part for = 6 or t is given by 

da^o^^igg ^ q'q') _ Trag + 3(^^ + t') 

dcosO As ^ 3ut 4s2 J ^ ' ^ 

and in a complete computation the SM QCD corrections will have to be added. 

As one sees from eq.(3.6,3.7), as a "technical" consequence of the t, u channel diagrams 
dominance in the chosen configuration, the SUSY electroweak and QCD effects concur at 
their (expectedly accurate) one-loop level to produce an overall negative effect that can be 
sizable, particularly for large tan /3 values where it could reach a relative twenty percent 
size, as we shall show in detail in the following Section. In this sense, and within the special 
scenario of large cm. energies and "reasonably" light SUSY masses that we have fixed in 
this simplified preliminary analysis, the chances of detecting a virtual SUSY effect in heavy 
quark production appear thus to be more promising for an experimental situation such 
that the considered initial gluon-gluon state gives (via its t and u channel diagrams) the 
dominant contributions to the cross section. From the available [16] luminosity pictures, 
we deduce that this request indicates the LHC experiments. Therefore, from now on we 
shall concentrate our investigation on this special machine, keeping in mind that for a 
different scenario, whose analysis might be performed in a less simple way, the role of 
TEVATRON might be relevant, or dominant, as well. 

As a final comment to this Section, we would like to remark the fact that the overall 
relative "genuine" SUSY correction at one loop, in the chosen LHC scenario, is almost 
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identical with that which one one would find, for the same heavy quark pair production in 
the identical cm. energy and SUSY scenario, at a lepton collider (LC). Without writing 
additional explicit formulae, we can simply explain this statement with the observation 
that, in the relevant gluon-gluon initiated process, the surviving SUSY effect is of purely 
universal (i.e. process independent) kind and due to the heavy quark final state, both 
for the clcctrowcak and for the QCD SUSY components, the latter ones being entirely of 
vertex kind owing to the previously shown RG suppression. These contributions factorize 
in the same way and therefore produce the same relative SUSY effect for initial gluon- 
gluon and electron-positron state, even if the relative Standard Model corrections may be 
different for the two cases e.g. when non universal angular dependent terms appear. The 
only (small) differences are due to the universal (non Yukawa) SUSY contribution (i.e. 
the — /n(s/M^) term) arising from the initial e+e" state. 

Our analysis of the simple partonic processes is thus completed. The next step is now 
that of considering to which amount the features that we have underlined will survive 
in the real process of production from a proton-proton state. This will be done in the 
forthcoming Section. 



IV. CROSS SECTION FOR HEAVY QUARK PAIR PRODUCTION IN PP 

COLLISIONS 



We now consider proton-proton collisions with inclusive production of a pair of heavy 

quarks PP — > q'q' + In this preliminary analysis we just want to show the role of 

the SUSY corrections on both quark- ant iquark and gluon-gluon subproccsscs. With this 
purpose we will concentrate our attention on the invariant mass distributions of final bb 
or tt quarks. Future works may consider other types of distributions (like Pt distributions 
of the quarks or of their decay products) using the subprocess amplitudes that we have 
established in Sect.II and III, and the corresponding parton model kinematical tools. 

For a total cm. squared energy S, the q'q' squared invariant mass (s) distribution is 
given by 

where t = ^, and (ij) represent all initial qq pairs with q — u,d, s, c, b and the initial gg 
pairs, with the corresponding luminosities 

L,jiT,cose) = — — /_ dy [ ^{x)J{-)+J{x)^{-) ] (4.2) 

i{x) being the distributions of the parton i inside the proton with a momentum fraction, 
X = e^, related to the rapidity y of the q'q' system. The limits of integrations for y 
can be written 

Vmax = max{0, min{y - ^Inx, Y + ^Inx, - ln{^/T)}} 

Umin — Umax 

(4.3) 
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where the maximal rapidity is y = 2, the quantity x is related to the scattering angle in 
the g'g' cm. 

1 + cos6' ,^ ^, 

X = -n 4.4 

1 — cos 6 

and 



,2 



COS dmin.max = T V 1 ^ (4-5) 

expressed in terms of the chosen value for PT,min- 

We have evaluated numerically the above expression of the differential cross section 
in the LHC case with = 14 TeV, a fixed SUSY mass scale Msusy = 350 GeV and 
an angular cut corresponding to pT,min = 10 GeV. Concerning the parton distributions, 
we must stress that in principle we should consider their evolution up to the desired 
energy scale y/s in the framework of SUSY QCD. However, the supersymmetric corrections 
to the evolution should be negligible in a first approximate treatment if the masses of 
the supersymmetric particles are large enough, still not spoiling the applicability of the 
Sudakov expansion. With these remarks in mind, we have used the 2003 NNLO MRST 
set of evolved parton distribution functions available on [16]. 

A summary of our numerical analysis is shown in Figs. (6-8). In Fig. (6), we show 
the percentual effect on the cross sections for production of final bottom or top quarks 
at two representative values tan/? = 10,40 in the cm. energy range 0.7-1 TeV (the lower 
limit corresponds, qualitatively, to a value s = 4M|^^y where we can still hope from our 
experience in H^H~ study [11] that our logarithmic expansion is "reasonable"). For the 
higher value, the effect reaches the remarkable 20 % level. An analysis of the relative 
weights of the various subprocesses contributing the total cross section shows that the 
dominant subprocess is the gg one. This is due to the low values of the ratio s/S dX LHC 
in the considered range for the final state invariant mass. Indeed, at low s/ S the fraction x 
is also typically small and the rapid rise of the gluon distribution function overwhelms the 
role of the other subprocesses; see for example the illustrations for gg and qq luminosities 
given in ref. [12]. To give some numerical examples one can check that for ^/s — 1 TeV 
the contribution of the gg subprocess represents about 80% (88%) of the total Born cross 
section for final bottom (top). At the smaller ^/s = 700 GeV, the effect is even larger 
with the gg subprocess being now 85% (92%) for final bottom (top). Since the process is 
dominated by the gg subprocess, the features of Figs. (6-8) can be explained in terms of 
Eqs. (3.6,3.7). In particular at the special value tan P — AO the two Yukawa combinations 
proportional to and ml appearing in Eqs. (3.6,3.7) happen to be equal explaining the 
almost superposed lines in the plot ^. 



We used in this preliminary analysis the fixed values mt = 173.8 GeV, = 4.25 GeV, 
rather than (energy dependent) running values. In so doing, we ignored a higher order effect 
consistently with the philosophy of our paper. 
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In Fig. (7), wc show the comparison between the effects in the Standard Model and 
those that we find in the MSSM at tan/i = 40, a vahie that has the "advantage" of 
providing a large correction due to the Yukawa terms. Indeed, it is precisely this kind 
of contribution that is responsible for the significant enhancement of the effect compared 
with the Standard Model case. The reason is the amplification of the coefficient of the 
term by more than a factor 2 due to the replacement — >■ 2mj(l + cot^/3) as 
well as the additional large correction tan^/3 introduced by the analogous replacement 
ml — > 2m^(l + tan^ j3). As a comment about the numerics we remark that in the Standard 
Model the Yukawa effect for final top is larger than in the case of final bottom by the 
factor (Sm^ + ml) / (rn^ + ml) ~ 3 explaining why the Standard Model full effect is larger 
for final top. In the MSSM at tan /? ~ 40 we already discussed the equivalence of the 
Yukawa effect for the final top or bottom. 

Finally, in Fig. (8), we show the relative weights of the genuine SUSY contributions 
that are not present in the Standard Model. They are three, i. e. (i) the SUSY component 
of the QCD correction (we already mentioned the genuine Standard Model QCD correction 
and we simply remark that is shared by both the Standard Model and the MSSM), 
(ii) the SUSY gauge electroweak Sudakov (linear) logarithmic term, (iii) the additional 
tan/3 dependent Yukawa terms. The dominant effects are the Yukawa and the SUSY 
QCD correction, the first one being, as anticipated in the Introduction, the largest. In 
agreement with the previous Figure (7), the Yukawa part of the difference "MSSM - SM" 
is larger for final bottom (at tan/9 = 40). Indeed, as we remarked, we have MSSM(top) 
~ MSSM(bottom) and SM(top) > SM(bottom) leading to MSSM-SM(top) < MSSM- 
SM(bottom). Notice also that minor numerical differences between the curves for final 
bottom or top must be traced back to the fact that for final bottom the cross section 
contains a small, but non negligible, component from the subprocess bb — > bb whose Born 
angular dependence is totally different than the counterpart bb — > tt in the case of final 
top. 

Figures (6-8) show the main result of our paper. One sees that the relative overall 
SUSY effect could be large, varying from approximately ten percent to approximately 
twenty percent in the range tan P — 10 — 40. This effect is definitely larger than the corre- 
sponding SM one. as shown by Fig. (7). In particular, one notices that the enhancement 
is less due to the SUSY QCD contribution, and is mostly coming from the SUSY Yukawa 
term, which would be absent in the case of light quark production. For an experimen- 
tal and theoretical precision at the few percent level (i.e. summing statistics, detection 
efficiencies and uncertainties in quark and gluon distributions), the presence of such a 
SUSY correction represents therefore, in our opinion, a feature of the process that cannot 
be ignored, and could provide a rather stringent test of the Supersymmetric model to be 
investigated. In this sense, the production of heavy quark pairs at LHC appears to us to 
be particularly interesting. 
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V. CONCLUSIONS 



A number of realistic statements must be made when drawing some possible conclu- 
sions from this paper. Our analysis has undoubtedly been specific, since it has assumed 
a combination of events i.e. a previous discovery of Supersymmetric particles and a "rea- 
sonably light " nature of the SUSY scenario. In this particular case, we have shown that 
in the production of heavy quark pairs at cm. energies in the one TcV range, where quite 
reasonably a one-loop logarithmic Sudakov expansion should provide an accurate descrip- 
tion of the virtual SUSY MSSM correction, the latter could reach values in the twenty 
percent range for large tan (5 and might be therefore detectable at proton colliders, in 
particular at the LHC which seems to us to be, for the specific cm. energy configuration 
that we have chosen, the more suitable machine. Of the relevant logarithmic expansion 
we have computed the leading (quadratic) and next to leading (linear)term, the SUSY 
genuine effect being simply of linear kind. We have performed this computation leaving 
undetermined a possible next-to-next-to leading term, in practice a constant one. In a 
complete treatment, the latter term should be computed or at least estimated. This is 
not trivial since in this quantity a large number of parameters of the supersymmetric 
model will generally appear, and several extra assumptions should be made concerning 
their values that are, we believe, beyond the realistic purposes of this first preliminary 
analysis ^. With this premise, it seems to us that a relevant feature that emerges is that 
the relative size of the genuine SUSY correction could be large, definitely larger that in 
the SM case, mostly owing to the role of the Yukawa correction. 

Another comment is related to the possibility that the SUSY scenario is not as "rea- 
sonably" light as we assumed. In the case for which the relevant masses are not huge, 
we feel, though, that the SUSY virtual effect, to be computed in a less simple way i.e. 
without logarithmic expansions, should still be numerically similar to the one that we 
computed, i.e. should not depend dramatically on the values of the SUSY masses of the 
process and should still be computable in a one loop approximation, and we will devote 
a future investigation to this special (negative) situation. 

For what concerns a comparison with other similar work, we remark that our conclu- 
sions concerning a large virtual effect in the MSSM arc in line (but with SUSY enhanced 
contributions) with those of an analysis of WW production in the Standard Model frame- 
work [18]. It should also be recalled that, again within the Standard Model framework, 
the production of b pairs appears to be promising for the detection of virtual effects in 
precision measurements [4]. In this sense, we believe to have shown that this conclusion 
could still be valid (and even more spectacular) for a more general class of final pairs, in 
a (hopefully valid) Supersymmetric picture. 



^In particular, given the fact that all the logarithmic SUSY effects at one loop are of Unear order, 
possible SUSY masses Mi larger than MsusY (e.g. masses of heavy gluinos) are automatically 
reabsorbed by constant terms ~ alogMi/MsuSY that should though, remain tolerably small, 
given the absence of "Yukawa enhancement" . 
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(d) (e) (f) 

FIG. 2. Diagrams for electroweak corrections to the annihilation amplitude qq q'q'; Stan- 
dard Model (a,b,c), where solid lines represent quarks, in (a,b) dashed lines represent gauge or 
Higgs bosons and in (c) the wavy line is a photon or a Z; SUSY diagrams (d,e,f), where in (d,e) 
the internal solid line is a gaugino and the dashed line is a squark, and in (f) solid lines represent 
quarks and the dashed line is a SUSY Higgs boson. 
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(c) (d) 




(e) 

FIG. 3. Diagrams for SUSY QCD corrections to the qq q'q' annihilation amphtude 
(a,b,c,d), in which the internal solid lines are gluinos and the dashed lines are squarks, and 
diagram for gluon self-energy diagrams (e). 
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(a) 

FIG. 4. Born diagrams for gg 
exchanges (b). 



(b) 



q'q' , s-channel gluon exchange(a) and t-,u- channel quark 
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(c) (d) 

FIG. 5. Diagrams for SUSY QCD corrections to gg — > q'q' in the t-,u- channels (a,b) (one 
should add similar diagrams to (a,b) with down triangles) and in s-channel (c,d); in (a,c) the 
triangles contains squark (dashed) and gluino (solid) lines; in (b,d) they contains quarks (solid) 
and SUSY Higgs bosons (dashed) lines. 
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FIG. 6. Effect of the combined electroweak and SUSY QCD corrections in the cross section 
for final bottom or top pairs at LHC. The various parameters are \fS = 14 TeV, MguSY = 
350 GeV and pT,min = 10 GeV. We show the results obtained with two values of the MSSM 
parameters tan/3. 
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FIG. 7. Comparison of the full effect in the Standard Model and in the MSSM. We show the 
correction to the cross section for final bottom or top pairs at LHC in the Standard Model and 
in the MSSM with tan (3 = 40. The other parameters are as in previous Figure: s/S = 14 TeV, 
MsuSY = 350 GeV and pT,min = 10 GeV. 
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FIG. 8. Separate SUSY effects. We show the three extra SUSY effects that arc present in 
the MSSM with respect to the Standard Model. They are the SUSY component of the QCD 
correction, the universal SUSY electroweak terms and the tan/3 dependent SUSY Yukawa terms 
computed at tan /? = 40. We recall that these genuine SUSY contributions grow like log s. The 
other parameters are a/5 = 14 TeV, MsuSY = 350 GeV and PT,min = 10 GeV. 
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